Control via electron count of the competition between magnetism and 
superconductivity in cobalt and nickel doped NaFeAs 
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Using a combination of neutron, muon and synchrotron techniques we show how the magnetic 
state in NaFeAs can be tuned into superconductivity by replacing Fe by either Co or Ni. Electron 
count is the dominant factor, since Ni-doping has double the effect of Co-doping for the same doping 
level. We follow the structural, magnetic and superconducting properties as a function of doping to 
show how the superconducting state evolves, concluding that the addition of 0.1 electrons per Fe 
atom is sufficient to traverse the superconducting domain, and that magnetic order coexists with 
superconductivity at doping levels less than 0.025 electrons per Fe atom. 



Since the discovery of high-temperature superconduc- 
tivity in fluoride-doped LaFeAsO p], several series of 
compounds containing fluorite-type or anti-PbO type 
iron pnictide or chalcogenide layers with a range of in- 
tervening layers have been shown to exhibit unconven- 
tional superconductivity [2J. The superconducting pnic- 
tides range from compounds in which the iron pnictide 
layers are separated by electropositive metal cations to 
those in which a much thicker perovskite-type oxide slab 
separates the iron pnictide layers and the electronic di- 
mensionality of the system is correspondingly much lower 
[3]. As for the layered cuprate superconductors, the su- 
perconducting regime is adjacent to antiferromagnetic or- 
der. The boundary between these regions may be tra- 
versed by changing the electron count, by applying an 
external hydrostatic pressure or by making an isoelec- 
tronic chemical substitution which changes the structural 
parameters, essentially applying a "chemical pressure". 

The chemically simplest superconductors containing 
iron arsenide layers have formula AFeAs where A=Li, 
Na and adopt the anti-PbFCl structure type [see the in- 
set to Fig. [TJb)]. LiFeAs was shown to be a bulk su- 
perconductor in its undoped, stoichiometric form [U [5] . 
Stoichiometric NaFeAs [6] also showed evidence for su- 
perconductivity: a sample which appeared to consist of 
a single crystallographic phase with very close to full oc- 
cupancy of the Na site, and which had a 10% supercon- 
ducting volume fraction evident from DC magnetometry, 
also exhibited long-range antiferromagnetic ordering in 
most of its volume as demonstrated by muon-spin rota- 
tion (/iSR) measurements [BJ. These results imply that 
there is some overlap between the antiferromagnctically 
ordered and superconducting regions of the phase dia- 
gram. In this Letter we present high-resolution powder 
diffraction, ^tSR, and magnetometry measurements on 
NaFei-^M^As (M = Co, Ni) that show the extremely 
high sensitivity of the competition between magnetic or- 



dering and superconductivity in NaFeAs derivatives to 
the electron count. 

The compounds prepared were NaFei-^Co^As [x = 
0, 0.005, 0.01, 0.025, 0.05, 0.07, 0.1, 0.15, 0.2) and 
NaFei-^NirAs [x = 0.025, 0.05, 0.1, 0.2). All are ex- 
tremely air sensitive and all manipulations of solids were 
performed in an argon-filled glove box. The products 
were prepared from stoichiometric quantities of elemen- 
tal reagents: freshly cut sodium pieces were placed in a 
9 mm diameter Nb tube and a well-ground mixture of 
the transition metal and arsenic powders were placed on 
top. The tube was sealed under 1 atm of argon in an 
arc- welding furnace and heated at 200° C in a resistance 
furnace for 48 hours. The black initial product was ex- 
tracted from the tube, homogenized in an agate mortar 
and pressed into a pellet which was heated in a second 
sealed Nb tube at 750° C for 48 hours. The tube was 
cooled to room temperature in a few hours at the natu- 
ral rate of the furnace. The final products were highly 
crystalline lustrous black powders and appeared phase 
pure according to the results of laboratory X-ray powder 
diffraction measurements. 

Structural investigations were carried out using the 
high-resolution X-ray diffractometer ID31 at the ESRF, 
Grenoble, France and the high-resolution neutron pow- 
der diffractometer HRPD at the ISIS facility, UK. On 
ID31 the samples were contained in 0.7 mm diameter 
glass capillary tubes sealed closed under helium exchange 
gas. Samples for neutron diffraction were contained in 
vanadium cylinders sealed with indium gaskets. Rietveld 
refinement was performed using the GSAS suite [7] . The 
diffraction data did not indicate more than 1-2% Na de- 
ficiency, consistent with earlier reports [SI [5] . The evolu- 
tion of the lattice parameters and key structural param- 
eters with composition at room temperature are summa- 
rized in Fig. [I] The Fe(Co)As4 tetrahedra become more 
regular as Fe is substituted by Co [Fig. [lja)] and the 
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FIG. 1: (Color online.) The evolution at 298 K of (a) As-Fe- 
As bond angles (b) lattice parameters and (c) Fe-As bond 
lengths of NaFei-zCo^As as a function of Co content. The 
lines are straight-line fits. Inset: crystal structure of NaFeAs. 



Fe(Co)-As bond distances decrease almost monotonically 
[Fig. [ljb)] as the smaller Co atom is substituted. The 
almost monotonic contraction in cell volume (by ~ 1% 
between x = and x = 0.2) is dominated by the contrac- 
tion of the c lattice parameter (« 0.8% over this range of 
compositions, compared to w 0.1% for the basal lattice 
parameter, a [Fig. |ljc)] ). A similar effect is found for Ni 
doping. 

Magnetometry measurements were carried out on a 
Quantum Design MPMS-XL SQUID magnetometer un- 
der zero-field cooled and field cooled conditions in a 
measuring field of 5 mT. The evolution of the low tem- 
perature magnetic susceptibility with composition is 
shown in Fig. [2j The reported behavior of NaFeAs 
13 EJ [10] is reproducible with several samples showing 
very small superconducting fractions (5-10%) and broad 
superconducting transitions with T c of around 10 K as 
judged from the onset of diamagnetism. Small levels of 
electron doping away from NaFeAs produce a rapidly- 
increasing superconducting fraction and in the case of 
Co doping, which supplies one additional electron per Co 
atom, a fully superconducting state is attained between 
NaFeo.975Coo.o25As and NaFeo.93Coo.07 As with the max- 
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FIG. 2: (Color online.) The evolution of the zero-field mag- 
netic susceptibility (measured in an applied magnetic field of 
5mT) of Nai-^CoxAs samples (left) and the comparison of 
Ni-and Co-doped samples (right) showing that the suscepti- 
bilities are similar for isoelectronic compounds. Doping levels 
are expressed as percentages of Co unless specified. 



imum T c of 21 K attained in NaFeo.97sCoo.o25As. 

Substituting Ni for Co at an equivalent doping level 
doubles the amount of electron doping but produces 
a very similar doping dependence of lattice param- 
eters. As shown in Fig. [2] the isoelectronic com- 
pounds NaFeo.95Coo.05As and NaFeo.975Nio.o25As ex- 
hibit extremely similar susceptibilities. Similarly, at the 
electron-rich boundary of the superconducting region, 
the isoelectronic NaFeo.9Coo.1As and NaFeo.95Nio.05As 
both show very tiny superconducting fractions. These 
results strongly suggest that electron count, and hence 
band filling, is the key parameter which controls the su- 
perconductivity in this series of compounds. Samples 
with higher doping levels than 0.1 electrons per Fe did 
not exhibit superconductivity. 

Zero-field /iSR experiments were carried out on the 
GPS beamline at the Swiss Muon Source (S^S). In 
NaFeAs, a clear oscillatory signal is observed [Fig. |3ja)] 
which can be fitted to an expression which is the sum 
of three components: A(t) = Y^i=i Ai cos(27Wj£)e _Ait . 
All three precession frequencies Vi can be fitted and each 
follows a temperature dependence [Fig. [3j^b)] consistent 
with a Neel temperature Tn = 45.0(2) K and a critical ex- 
ponent f3 = 0.20(2); These results are consistent with an- 
tiferromagnetic order with a primarily two-dimensional 
character. The relaxation rates Aj stay roughly constant 
with temperature, but one of them (corresponding to the 
highest frequency component) diverges at Tn. On substi- 
tuting 1% of the Fe ions with Co to obtain a sample with 
a superconducting fraction of about 60%, as judged by 
magnetometry, and a broad superconducting transition 
[Fig. [2] , there is still evidence for long-range magnetic or- 
dering in the sample, but the oscillations in the asymme- 
try are heavily damped [Fig.[3]jd)] and lower in frequency 
and Tn [Fig. kjje)]. This implies a smaller internal field 
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FIG. 3: (Color online.) (a) Zero-field /xSR data in NaFeAs 
at 1.8 K. (b) Temperature dependence of the three preces- 
sion frequencies (data for two different samples are shown in 
open and closed symbols), (c) Relaxation rate for the three 
frequencies, (d) Zero-field piSR data in NaFeo.99Coo.01 As at 
1.8 K and (e) temperature dependence of the precession fre- 
quencies (only two frequencies can be extracted in this case). 



with a more inhomogeneous field distribution. 

The undoped and lightly Co-doped samples were mea- 
sured using high-resolution diffraction measurements at 
low temperatures. The non-superconducting "parent" 
phases of these arsenide superconductors commonly show 
a structural distortion driven by the long-range anti- 
ferromagnetic ordering. Such a distortion was not ev- 
ident in the original report on NaFeAs [6], but it was 
observed in a higher resolution study [5] and the sym- 
metry of the distortion was observed to be similar to 
that in the "1111" and "122" parent phases with an or- 
thorhombic v2flT X V^ot X ct cell of Cmma symmetry 
adopted at low temperatures where g,t and ct are the 
tetragonal lattice parameters. On HRPD the splitting 
of, for example, the 112 reflection in P4/nmm into the 
022 and 202 in Cmma is well resolved [Fig. |4ja)-(c)]. 
The high-resolution measurements show that in NaFeAs 
the distortion does not vanish until 55-60 K, well above 
the antiferromagnetic Neel temperature similar to the be- 
haviour observed for all classes of parent materials apart 
from the AFe 2 As 2 (v4=Ca, Sr, Ba) "122" series where 
the two transitions are coincident. Our measurements 
confirm that the two transitions identified by Chen et al. 
[9] at 41 K and 52 K from heat capacity measurements 
correspond to the antiferromagnetic-ordering transition 
at 45 K and the structural distortion at about 55 K. 

In the lightly Co-doped regime x < 0.01 in which ^SR 
shows that antiferromagnetic ordering persists, diffrac- 
tion measurements on the high-resolution powder diffrac- 
tometer ID31 show that the distortion (bo/ao — 1, in 
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FIG. 4: (Color online.) (a) Rietveld refinements of NaFeAs 
against 2K neutron data (HRPD). Data (red, crosses), fit 
(green solid line) and difference plots are shown. Scattering 
from the sample environment has been excluded, (b)-(c) De- 
tail showing the splitting of the tetragonal 112 reflection into 
the orthorhombic 022 and 202 reflections, (d) Evolution of 
the orthorhombic 022/202 reflections in NaFeo.ggsCoo.oosAs 
measured on ID31; At 46 K the orthorhombic and tetrago- 
nal models give fits of comparable quality, (e) Temperature 
dependence of the distortion for NaFei-zCo^As. The error 
bars lie within the points except when the presence of the 
distortion is marginal, (f ) The size of the distortion at 5 K as 
a function of composition. Inset: for NaFeo.97sCoo.o25As the 
tetragonal 212 reflection measured on ID31 at 5 K is broader 
than at room temperature while the tetragonal 104 reflection 
(orthorhombic 114) does not broaden. 



terms of the basal plane orthorhombic lattice parame- 
ters) persists above T N [Fig. Hfd)-(e)]. The size of the 
distortion at 5 K decreases smoothly with increasing Co 
content [Fig. |4^f )] showing a similar dependence on com- 
position to the Neel temperature. 

In iron pnictide systems the structural distortion as- 
sociated with antiferromagnetic order is known to be 
quenched in the fully superconducting regime. Measure- 
ments at 5 K on ID31 of NaFeo.975Coo.025 As, which lies in 
this regime, revealed broadening of the reflections which 
would split under an orthorhombic distortion, while re- 
flections which would not split retained their ambient 
temperature widths [Fig. |4^f) inset]. At 5K, the or- 
thorhombic model refined stably but the goodness of fit 
was negligibly superior to that of the tetragonal model 
and the lattice parameters were clearly correlated with 
the peak profile parameters. An orthorhombic distortion 
may thus persist into the superconducting regime, al- 
though the upper bound on bo /oq — 1 is 3.0(1) x 10~ 4 in 
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FIG. 5: (Color online.) Uemura plot of the superconducting 
transition temperature T c versus the low temperature super- 
fluid stiffness. Data obtained here for the doped NaFeAs sam- 
ples are compared with previously reported data for other iron 
arsenides and cuprates [12H15] . The lower dashed line indi- 
cates the trend line for electron-doped cuprates, which closely 
corresponds to the behavior of the LiFeAs system [16] , in con- 
trast to the other pnictide superconductors, which sit closer 
to the trend for hole-doped cuprates. The symbol size rep- 
resents the typical estimated uncertainty of the data points. 



Ps using Ps /c = Xj = (3/0.0Q37l)V 2 B IIBS /<f>o- Plotting 
our data on an Uemura plot (Fig. [5]) demonstrates that 
the values of superfluid stiffness correlate well with the 
trend observed in hole-doped cuprates and other iron ar- 
senides, in contrast to the enhanced superfluid stiffness 
found for LiFeAs [IB]. 

The conclusions of our work are that NaFeAs lies very 
close to one edge of a superconducting dome which is 
traversed by the addition of 0.1 electrons per Fe atom 
using either Co or Ni doping (see Fig. [6]). As the sys- 
tem passes from the magnetically ordered regime to the 
superconducting regime there is a region of phase coexis- 
tence which has a width of less than 0.025 electrons per 
Fe, so the system is extremely sensitive to composition. 
The width and shape of the superconducting dome and 
the width of the coexistence region are similar to those 
determined in the "122" system Ba(Fei_ a; Co x )2As2 [17] 
and in the "1111" system CaFe 1 _ x Co x AsF [IB] . 
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assistance on HRPD, and A. N. Fitch and I. Margiolaki 
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at the Swiss Muon Source (SpS). 
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FIG. 6: (Color online.) Phase diagram for NaFei-^Coa; As. 
The maximum temperature at which the distortion persists 
is shown by the filled circles and dotted line. The antiferro- 
magnetic (AF) phase is delineated by Tn (open circles) and 
superconducting (SC) phase by T c (filled squares). 



NaFeo.975Coo.o25As, an order of magnitude smaller than 
in NaFeAs, and marginal even using an extremely high- 
resolution diffractometer. 

Once a large superconducting fraction is established 
(x > 0.025) it is possible to use transverse-field /iSR to 
study the superfluid stiffness by measuring the temper- 
ature and field-dependence of the rms field broadening 
£> rms due to the vortex lattice |llj . Such experiments 
were performed using both S/xS and the ISIS Pulsed 
Muon Facility and these data can be used to deduce the 
in-plane penetration depth A Q {, and superfluid stiffness 
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